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ABSTRACT 

A systematic investigation on the evolution of a prolate cloud at an Hu boundary is 
conducted using Smoothed Particle Hydrodynamics (SPH) in order to understand the 
mechanism for a variety of irregular morphological structures found at the boundaries 
of various Hu regions. The prolate molecular clouds in this investigation are set with 
their semi-major axes at inclinations between 0 and 90° to a plane parallel ionizing 
radiation flux. A set of 4 parameters, the number density n, the ratio of major to minor 
axis 7 , the inclination angle ip and the incident flux Peuv, are used to define the initial 
state of the simulated clouds. The dependence of the evolution of a prolate cloud under 
Radiation Driven Implosion (RDI) on each of the four parameters is investigated. It is 
found that: i) in addition to the well studied standard type A, B or C Bright Rimmed 
Clouds (BRCs), many other types such as asymmetrical BRCs, filamentary structures 
and irregular horse-head structures could also be developed at Hu boundaries with 
only simple initial conditions; ii) the final morphological structures are very sensitive 
to the 4 initial parameters, especially to the initial density and the inclination; iii) 
The previously defined ionizing radiation penetration depth can still be used as a 
good indicator of the final morphology. 

Based on the simulation results, the efficiency of the RDI triggered star formation 
from clouds of different initial conditions is also estimated. Finally a unified mechanism 
for the various morphological structures found in many different Hu boundaries is 
suggested. 

Key words: Hll regions - hydrodynamics - stars: formation - ISM: evolution - ISM: 
kinematics and dynamics - radiative transfer 


1 INTRODUCTION 

As the intensive EUV radiation from newly formed star(s) 
ionizes the star-facing surface layer of a molecular cloud, the 
ionization heating induces a shock which propagates into 
the cloud, compressing material to form highly condensed 
core(s). These cores are the potential sites for EUV radia¬ 
tion triggered new star formation. Simultaneously the Ha 
emission produced by the recombination of ions with elec¬ 
trons creates a bright rim around the star facing side of the 
cloud. The structures formed are termed Bright Rimmed 
Clouds (BRCs), the best candidate for studying the feed¬ 
back of massive stars to surrounding and parental molecu¬ 
lar cloud (s). This process is the Radiation Driven Implosion 
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(RDI) mechanism for EUV radiation triggered star forma¬ 
tion in BRCs. 


Multi-wavelength observations have revealed various 
physical and morphological properties of BRCs at Hll re¬ 
gions. BRCs observed so far tend to be classified as type 
A, B and C in increasing order of the curvature of their 
bright rim (Sugitani et al. 1991| Sugitani & Ogura 1994), 
though some of BRCs are found in M shaped morphology 


( Osterbrock||1957 Urquhart et al.||2006 Karr et al.|2005 ). 


Starting with a uniform and spherical molecular cloud 
illuminated by a plane-parallel ionizing radiation from one 
side, theoretical modelling based on the RDI mechanism suc¬ 
cessfully revealed the formation process of standard type A, 
B or C BRCs. These simulated BRCs are all symmetric to 
their structure axis, and the latter aligns with the radia¬ 
tion flux direction, i.e., the line connecting the tip of a BRC 
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Figure 1 . The distribution of the inclination angle which is 
defined as the angle between the major axis and the x axis as 
shown in Figure |2| of the BRCs in about 10 emission nebulae 
jOsterbrock|1957| K 


and the centre of the exciting star (|Bertoldi 1989| |Lefloch 

& Lazareff 

1994; Kessel-Deynet & Burkert 2003;; Miao et al. 

2006| 

2009 

Gritschneder et al.I2009;|Miao et al. 2010; Bisbas 

et al.| 

2011 

Haworth & Harries|2012| Tremblin et al.|2012). 


However, further investigation on the observed struc¬ 
tures at Hu boundaries finds that: i) some BRCs do not show 
a comet-like morphology, but fragment-core linear struc¬ 
tures perpendicular to the EUV radiation flux direction (e.g. 


Chauhan et al. (2011)); ii) most type A, B and C BRCs 


are not symmetric about their structural axes, which also 
do not necessarily align with the radiation flux direction of 


their exciting stars (Thompson et al. 2004 Morgan et al. 
2004). As shown in Figure [lj there is a distribution in the 


angle between the axis of a BRC and the direction from the 
tip of the BRC to the centre of the exciting star, based on 
the observational data on groups of BRCs in about 10 emis¬ 
sion nebulae (Osterb rock||1957 ); iii) There are more objects 
with structures which can not be categorized by type A, B 
or C than those which can be; iv) In those asymmetrical 
BRC structures, RDI triggered multi-star or sequential for¬ 
mation are often found not only at the head of BRCs, but 
also in the more compre ssed side layer (|Fukuda et al.||2013| 


Makela fe Haikala||2013| |Sicilia- Aguilar et al.||2014| |P anwar 

et al.|2014 ). It is obvious that it is difficult to accommodate 

the variety of structural and physical features observed with 
RDI modelling starting with a spherical cloud. 

Observations on the physical properties of molecular 
clouds have revealed that spherical clouds are very rare cases 


(Gammie et al. 2003 Doty et al. 2005 Rathborne et al. 2009) 


and there are many physical processes which lead to for¬ 
mation of prolate molecular clouds ( Gong k, Qstriker~||2011 


Gholipour & Nejad-Asghar 2013). To investigate the mecha¬ 
nism for the development of structures of various morpholo¬ 
gies other than standard (i.e., symmetrical) BRCs at Hu 
boundaries, using an initially non-spherical cloud, such as 
a prolate cloud, in RDI modelling seems the most feasible 
choice. Definition in terms of two additional, geometric, pa¬ 
rameters are introduced: the ratio of the major (a) to minor 
( b ) axis 7 = f and the inclination angle ip, the configuration 
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Figure 2. The initial geometry of a prolate cloud in simulations. 
x p , x a and x max are the x projections of the apex, the semi-major 
axis a and the farthest point of the cloud from z axis respectively. 


of which is illustrated in Figure [2] In this way we can ex¬ 
pect that the variety of the morphological structures derived 
from RDI simulations must be greatly increased. 

A set of RDI simulations were conducted using an ini¬ 
tially prolate cloud with p — 0°, which we call ‘perpen¬ 
dicular’ to the EUV radiation flux. It is found that a per¬ 
pendicular prolate cloud would evolve into a perpendicular 
fragment-core linear structure to the ionizing flux, covered 
by a bright top layer on the star side. Both the initial phys¬ 
ical/geometrical conditions of the cloud and the strength 
of radiation flux affect the distribution of condensed cores 
over the structure, the total condensed core mass and the 
triggered star formation time. The sporadic fragment-core 
structures at Hu boundaries can be interpreted by the ef¬ 
fect of the interaction between a perpendicular prolate cloud 
and a parallel-plane EUV radiation field. More details are 
presented in Kinnear et al. (2014). 

In this paper, we expand the investigation to the evo¬ 
lution of a prolate cloud at Hu boundaries, with the cloud’s 
major axis inclined to the EUV radiation flux by an angle 
0° ^ (90° — <£?)< 90°. This is performed in order to under¬ 
stand the mechanism for the formation of a broad range of 
asymmetrical structures and the consequent triggered star 
formation inside them. In the rest of the paper, we first 
present a brief description of the codes used, the initial con¬ 
ditions of the clouds for simulations, and a derived quantity 
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which gives an indication of the dynamic evolution of pro¬ 
late clouds at different Hu boundaries. In the results and 
discussion section we present simulation results and discuss 
the physical mechanism responsible for the formation of dif¬ 
ferent mophological structures at Hu boundaries. Finally, 
conclusions are drawn. 


2 THE CODES AND INITIAL CONDITIONS 


2.1 The codes 

An updated Smoothed Particles Hydrodynamics (SPH) 
Code ii is used for all of the simulations investigated in 
this paper, which is an extended version of Code I ( |Nelson 


V Langer 1999). Code i is an SPH code to investigate the 


effect of an isotropic FUV interstellar background radiation 
on the evolution of a molecular cloud, by solving the full 
set of hydrodynamic equations, including self-gravity of the 
cloud and a chain of chemical differential equations. Based 
on Code i, a numerical solver with a ray-tracing algorithm 


(Kessel-Deynet V Burkert 2000) for a plane-parallel EUV 


radiation transferring equation was implemented, as well as 
the relevant heating and cooling processes. We refer to this 
version as Code ii, which can be used to investigate the in¬ 
teraction of a plane-parallel ionizing radiation from massive 
star(s) with a molecular cloud of arbitrary initial geometry. 
More detailed descriptions of the code can be found in Miao 
|et al. ( |2006 ) and Kinnear et al. (2014). 

By virtue of the asymmetrical morphology of the simu¬ 
lated clouds, the condensed cores triggered by RDI induced 
shocks are not usually aligned with the structural axis and 
multi-star formation is also frequently observed. To identity 
the positions of these triggered condensed cores and anal¬ 
yse their physical properties, a code called ‘CoreFinder’ is 
used, for which a detailed description is included in|Kinnear| 
etaLlpOlU l. For the simulations presented in this paper, all 
Tores’ are defined as a region of at least 0.03 M 0 (> 100 
SPH particles) and the density of e ach SPH partic l e sam - 


pled must be n ^ 10 b cm . As in Kinnear et al. (2014), 


the cores formed from the simulations are the first gener¬ 
ation of RDI triggered objects. Simulations for the subse¬ 
quent generations of RDI triggered star formation becomes 
excessively slower after initial extremely high density core(s) 
form, which leads to a decrease towards infinitesimally small 
time steps, and the effective halt of the simulation. This can 
be avoided by the use of the ‘sink’ particles of Bate &; Burk-| 


ert 


(1997), but is not currently implemented in Code ii. 


All of the column density images are produced using the 
software SPLASH, which is specially designed for processing 
SPH numerical data ( Price|| 2007). 


2.2 Initial conditions 

2.2.1 Initial geometry 

The geometrical shape of a prolate cloud in the simulations 
is described by a pair of initial parameters ( a, 7 ) and its 
orientation to the plane-parallel ionization radiation flux is 
0° ^ p ^ 90°, as shown in Figure [ 5 ] The ionization flux 
onto the surface of the cloud on the star side is treated as 
plane-parallel to the z axis, for an assumed large ratio of 
distance to the illuminating star against cloud size. 


In Figure [2] x p is the x coordinate of the apex of the 
cloud, letting r be the distance of any a point on the surface 
of a prolate cloud in the xz plane, and 0 the angle between 
r and the x axis, we have, 


* — —, — . v-l; 

yj cos 2 (0 — ip) + y 2 sin 2 (0 — c p) 

The maximum x coordinate and the x projection of the semi¬ 
major axis 


sin 2 p 


£max — a\/ COS 2 p + ^ 

a cos p. 


( 2 ) 

(3) 


As an important indicator to the dynamical evolution 
of a prolate cloud with an inclination angle <p, the effective 
illuminated area by the EUV radiation is, 


A = 


TT^max 

7 


— -- 4 / cos 2 p + 

7 


sin 2 p 


(4) 


which is a decreasing function with <p, has the maximum 
7 vab at p = 0° and the minimum Tvb 2 at p — 90°. For the 
convenience of the following discussion, we call the half cloud 
with p ^ 0 ^ tv + p as the ‘front’ semi-ellipsoid, and tv + p < 
6 ^2 tv + p the ‘back’ semi-ellipsoid. 


2.2.2 Initial mass distribution 


As we have done in the investigation of the evolution of a 


perpendicular prolate cloud at an Hu boundary (Kinnear 
eTaTpO lIt, all of the molecular clouds in our simulations 
start with a uniform density, which is rendered by a glass¬ 
like distribution of SPH particles created using Gadget-2 
( Springel] 2005). A glass-like distribution is taken as a good 
approximation to a uniform mass distribution. We choose 
the initial mass of the cloud to be 30 M 0 , the same as used 


for Kinnear et al. (2014). The number of SPH particles for 


each molecular cloud is determined by the mass resolution, 
3.0 x 10 -4 M 0 per SPH particle, for 100,000 particles, higher 
than required by the convergence test of the Code ii. A zero 
initial velocity held is set for all of the molecular clouds in 
the simulations. 

Since the objective of our investigation is to explore the 
evolution of the RDI triggered collapse of an initially in¬ 
clined prolate cloud, any initially unstable cloud against the 
self-gravity before interacting with the ionization radiation 
flux should not be chosen. Stability is assured by applying 
the Jeans criteria for an isolated prolate cloud, in terms of 
Jeans number J, the ratio of the initial gravitational en¬ 
ergy of a prolate cloud to the thermal energy and can be 
expressed as (Bastien 1983), 

J= *? P J lh l In 


15 e R g T 


m) 


s: 1.0, 


(5) 


where p, T and p are the mass density, the initial temper¬ 
ature and the mean molecular mass of the prolate cloud 
respectively, G and R g are the standard physical constants 
and the eccentricity, 


( 6 ) 


a z 7 

Using the relation between density, mass M and volume of 
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a cloud, we can derive the condition for the major axis a of 
an initially stable prolate cloud, 


CL ^ Ucrit — 0.052 


M* 7 

T \[l 2 - 1 


In 


( 7 + yV ~ 1 \ 
\7 - \J 7 2 - 1 ) ’ 


( 7 ) 


where M* is the mass of the prolate cloud in units of solar 
masses, and a and a cr it have units of Parsecs. For a given 
molecular cloud of mass M*, and initial temperature T and 
7 , a minimum value a can be estimated, the major axis of an 
initially gravitationally stable cloud should satisfy a > a cr it • 


2.3 Boundary condition 


The clouds in our simulations are subject to an isotropic 
interstellar background FUV radiation of one Habing unit 
(Habing 1968]) and an ionizing EUV radiation with a flux of 
lCrcni s -1 (typical of the boundary of an Hu region) di¬ 


rected parallel to the z -axis (along the negative z direction) 
as illustrated in Figure [2] 

A constant pressure boundary condition is applied with 
the value of the external pressure being set equivalent to an 
external medium of atomic hydrogen of n(Hi) = 10 cm -3 
and temperature of 100 K. An outflow condition is imposed 
at a fixed boundary which is a few times of the initial size 
of the cloud ( Nelson &; Langer|1999 ). 


2.4 The ionization penetration depth parameter 


In order to classify the initial dynamic regime of a perpen¬ 
dicular prolate cloud {jp = 0 ), a dimensionless parameter, 
cfeuv, was defined, being the ratio between the physical ion¬ 
ising radiation penetration depth to the semi-minor axis of 


the cloud. For this scenario, cfeuv = Feuv 1 (Kinnear et al. 

5 a :a n z '_ 

2014). When p ^ 0°, we modify its definition to the ratio of 


the physical ionization penetration depth to the character¬ 
istic depth of the cloud. This characteristic depth is defined 
as half of the longest path through the cloud in the direction 
of the radiation, which is always the depth of the cloud at 
x — 0 , y = 0 . 


ofeuv — 


^EUV 

n 2 


r(0 — 90°) 


Feuv 7 
as an 2 


= 1.6x10 


sin 2 p 


,3 Feuv 7 


+ cos 2 p 


sin 2 p 


+ cos 2 p (pc) ( 8 ) 


where the major axis a* is in the unit of pc, and uv is 
the EUV ionising radiation flux in units of 10 9 cm _ 2 s _1 , «b 
is the recombination coefficient of hydrogen ion - election 


under the ‘on-the-spot’ approximation (Dyson & Williams 
1997) and has the value of 2.0 x 10 ~ 13 cnU s -1 at a tem- 


perature of about 10 4 K ( Dyson V Williams||l997 ). This is 
then taken as a constant, as the equilibrium temperature for 
ionized material is « 10 4 K and the dependence of «b on 
temperature is not strong in the region around that temper¬ 
ature. 


As already discussed in Kinnear et al. (2014), when 


Feuv << 1 for p — 0, the prolate cloud is in the RDI trig¬ 
gered shock dominant region, and the collapse of the cloud is 



n 

7 

<P 

F 0 

Series name 

cm -3 


0 

xlO 9 cm -2 s -1 

G 1200 

1200 

1-10 

0-90 

1 

density & flux 

100-1200 

2 

60 

1 & 2 

irregular 

400 

2-3.5 

60-85 

2 


Table 1. Summary of the test series conducted, n is initial den¬ 
sity, 7 the axial ratio, p the inclination angle and Fo the incident 
flux. 


by ‘foci convergence’ mode, i.e., two high density cores form 
at the two ends of a filament; when Fe uv ^ 1 , the cloud can 
still be RDI triggered to collapse but through ‘linear conver¬ 
gence’, i.e., a few high density cores form along the whole 
filament; when Feuv ^ 1 , the cloud is in photo-ionization 
dominant region, and would be photo-evaporated. As Feuv 
can possess a similar range of values for all (p in Equation [ 8 ] 
we would expect it to play a similar role for p ^ 0 . 


3 RESULTS AND DISCUSSIONS 


Three sets of simulations were conducted using prolate 
clouds of different geometrical shapes ( 7 ), inclination an¬ 
gles (p), initial densities (n) and ionising fluxes (Feuv), to 
investigate the effect of their variations on the evolution of 
a prolate cloud. Table [I] contains a summary of the ranges 
of properties for each set. All of the clouds have an initial 
temperature of 60 K and are illuminated by a 50,000 K star 
providing a flux at the cloud, F 0 = 10 9 cm -2 s -1 (or 2.0 
when specified), the same as that for 


Kinn¬ 


ear et al. (2014). In this section, we first describe the evolu¬ 


tion of a prolate cloud of an inclination angle of 45° in detail, 
then investigate the influence of changing initial inclination 
angle p on the evolution of the same cloud. Next we dis¬ 
cuss how the initial geometry of a prolate cloud affects the 
morphological evolution of a cloud and RDI triggered star 
formation inside it. Finally we propose a formation mecha¬ 
nism for the variety in morphological structures found at H 
11 boundaries. 


3.1 Evolution of cloud G1200(5) with p s= 45° 

Table [2] lists all of the relevant parameters of clouds with 
initial density 1200 cm -3 and inclination angle p = 45°, for 
variations of the initial ratio 7 . The cloud is the fifth one 
in the G1200 set, notated G1200(5). It has an initial shape 
defined by 7 = 2.0, its major axis (0.784 pc) is ~ 10 times of 
its critical major axis (0.079 pc) derived by using Equation[7] 
indicating its initial stability against gravitational collapse. 

Figure [3] shows an evolutionary sequence of column den¬ 
sity of the cloud over a period of 0.13 Myr. We also plot the 
corresponding mean axial density distribution along the x 
axis in Figure [4j with the established binning method for 
describing axial profiles of SPH parti cle properties (|Kinn- 
ear et al. |2014[ |Nelson fc Langer|1999[ [Nelson fc Papaloizou 

1993). This view provides a qualitative picture of the loca¬ 
tion of condensed cores formed in the shocked layer of gas 
in the cloud by the RDI effect. 

It is seen from the above two figures that 28 Kyr after 
the radiation flux was switched on, a shock is established by 
the ionization heating. This shock starts compressing the 
cloud surface on the star-facing side. A slightly compressed 








































Evolution of Prolate Clouds at Hu Boundaries 


5 


G1200 


No. 

7 

1 

1.00 

2 

1.25 

3 

1.50 

4 

1.75 

5 

2.00 

6 

2.25 

7 

2.50 

8 

2.75 

9 

3.00 

10 

3.25 

11 

3.50 

12 

3.75 

13 

4.00 

14 

4.50 

15 

5.00 

16 

5.50 

17 

6.00 

18 

7.00 

19 

8.00 


«crit (pc) 
0.052 
0.060 
0.067 
0.073 
0.079 
0.084 
0.039 
0.093 
0.097 
0.101 
0.104 
0.107 
0.110 
0.116 
0.121 
0.126 
0.130 
0.138 
0.145 


«1200 (pc) 
0.494 
0.573 
0.648 
0.718 
0.784 
0.849 
0.910 
0.970 
1.028 
1.084 
1.139 
1.193 
1.245 
1.347 
1.445 
1.540 
1.632 
1.808 
1.977 


4uv (%) 
0.225 
0.218 
0.219 
0.221 
0.224 
0.228 
0.232 
0.237 
0.241 
0.246 
0.250 
0.256 
0.260 
0.268 
0.277 
0.285 
0.292 
0.307 
0.320 


Table 2. Parameters of molecular clouds of G1200 series with 
an inclination angle ip = 45°. From left to right, columns 1-3 are 
the number identity, axial ratio and the critical semi major axis 
defined by Equation [TlColumns 4-5 are the major axis and c?euv 
defined by Equation[8] All of the semi-major axes and critical 
semi-major axes are in units of pc and the penetration depth is 
in %. 


thin layer of mean density « 3.0 x 10 3 cm -3 appears at the 
EUV exposed surface, as shown in the upper-middle panels 
in Figures [3] and [4] The density at x = 0.49 pc is slightly 
higher than the axial mean density, which reflects the very 
early stage of a condensed core formation around the apex 
of the cloud. 

From t = 0.028 to 0.11 Myr, the shock continues propa¬ 
gating into the cloud. The density in the thin shocked layer 
continues increasing and the structure of the head becomes 
more distinctive. The non-head mean axial density and the 
central density of the head increase to 10 4 and 10 5 cm -3 
respectively at 0.054 Myr, and to 2.0 x 10 4 and 10 6 cm -3 
respectively at 0.08 Myr. At t — 0.11 Myr, the density in the 
shocked layer increases dramatically towards the apex of the 
cloud (at x — 0.32 pc). This density jumps from « 8 x 10 4 to 
8 X 10 5 cm -3 at x = 0.24 pc, then to the maximum 4 x 10 6 
cm -3 at the apex. It is interesting to find from the magnified 
image in the left panel of Figure [5] that a filament structure 
forms at the head, aligning with the EUV radiation direc¬ 
tion. The structure has a length of roughly 0.03 pc and the 
top end of it coincides with the apex of the cloud. 

The gas contained in both the shocked surface layer and 
in the filament becomes further compressed by t — 0.13 Myr. 
At this stage the maximum mean density in the shocked 
layer and in the filament are 10 5 — 10 6 and > 10 9 cm -3 
respectively. The right panel of Figure [5] further reveals 
the formation of two high density cores in close proximity. 
These dual cores may be the result of fragmentation of the 
filament-like structure observed at t = 0.11 Myr. 

Overall, the cloud G1200(5) with an initial inclination 
angle of 45° evolves to an asymmetrical type B BRC, with 
double RDI triggered high density cores embedded at its 
head, very close to the apex of the BRC. The cores are the 
potential sites for new star formation. At the same time 


0.15 


0.1 


0.05 


0 


- 0.05 



24 


22 


20 


0.26 0.28 0.3 0.32 0.34 0.26 0.28 0.3 0.32 0.34 

x parsec x parsec 

Figure 5. Magnification of the column density profiles time t = 
0.11 (left) and 0.13 Myr (right) at the head of the cloud G1200(5) 
with p = 45° 


v(°) 

4uv 

(%) 


0 

0.283 


15 

0.276 


30 

0.255 


45 

0.224 


60 

0.187 


75 

0.155 


90 

0.141 


Table 3. Variation of c/euv over different p for cloud G1200(5). 


the apex of the cloud moves down from z — 0.61 pc to 
z = 0.035 pc, a net recoil velocity of (v z ) — 4.3 km s _1 , 
due to the rocket effect caused by the evaporating gas from 
the star-facing surface irradiated by the ionizing radiation 
(Oort V S pitzer| 1955). Next we investigate the effect of the 
inclination of a prolate cloud on its dynamical evolution. 


3.2 Evolution of G1200(5) with 0° ^ (p ^ 90° 

To investigate the dependence on the inclination, the cloud 
G1200(5) is simulated for a range of inclinations 0° ^ ip ^ 
90°. The sequence of inclination angles are chosen as 0, 15, 
30, 45, 60, 75 and 90 °. The values of cfeuv for these seven 
clouds are listed in Table [3j from which we can see that 
the ionizing penetration depth decreases with p, and is very 
shallow, with values << 1. As a result, ah of these clouds are 
in the RDI induced shock-dominated regime and triggered 
collapse is expected. 


3.2.1 Morphological variation with p 

It is found that the evolutionary sequences of the 7 clouds 
are qualitatively similar to that with p = 45°, i.e., the for¬ 
mation of a condensed layer on the star facing surface of the 
cloud and highly condensed core(s) near the apex, although 
the final morphologies themselves depend on the initial incli¬ 
nation angle of the cloud. As the general evolution is so sim¬ 
ilar, they are not discussed individually, but with focus on 
comparison of the differences in their final morphologies and 
the physical properties of any high density core(s). In Figure 
[6] the final snapshots of column density of cloud G1200(5) 
from 6 different simulations with p =s 0,15,30,60,75, and 
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Figure 3. Sequential evolution of the column number density for Cloud G1200(5) with the incline angle (p = 45°. Time sequence is left 
to right, then top to bottom then left to right. 



x-axis (pc) 


Figure 4. Axial mean density distribution along the x axis of the cloud G1200(5) corresponding to the column density evolutionary 
snapshots in Figure [3] 
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90°, are presented. The equivalent for tp — 45° can be found 
in Figure [3) 

When the inclination angle ip is small 30°), x a ~ 
iCmax, the EUV radiation mainly illuminates the surface of 
the front semi-ellipsoid. The effect of the RDI is to drive a 
pseudo-plane-parallel shock propagating into the cloud, and 
the prolate cloud evolves to a filamentary structure with 
a condensed surface layer on the star side and highly con¬ 
densed core(s) at the apex (when p ^ 0°) or at either the 
apex or the two ends (when p = 0°). For the latter case, 
the evolution has been analysed in detail by[Kinnear et al.| 
(2014). 


When 45° ^ p < 90°, the EUV radiation illuminates 
not only the surface of the front semi-ellipsoid, but also the 
top part of back semi-elliptical surface, so that the shock in¬ 
duced by RDI creates a condensed and curved surface layer 
with the apex as a convergence point of two compressed 
curving surface layers, those surfaces on opposing sides of 
the apex. The region around the apex is also the site for the 
formation of high density core(s). These clouds evolve to 
asymmetrical BRCs of different morphologies. The morpho¬ 
logical sequence is similar to that of the cloud with p — 45°, 
as shown in Figure [3] Comparing the three final morpholo¬ 
gies of the simulations at 45, 60 and 75° in Figures [3] and 
[6] we can see that as the inclination angle increases, the 
final structure of the asymmetrical BRC becomes increas¬ 
ingly extended both in the head and tail. This is because, as 
the effective EUV illumination area decreases as described 
in Equation [4j the total volume of the cloud affected by 
the RDI induced shock becomes smaller. Consequently the 
original prolate shape of a cloud is less distorted. 

When p = 90°, the cloud G1200(5) develops into a sym¬ 
metrical type B/C BRC but having a wide tail structure, 
with the widest part defined by the initial semi-minor axis. 
The morphology developed is different from that of standard 
type B/C BRC structure, which has a much narrower width 
in its tail, if one assumes that the progenitor cloud was ini¬ 
tially spherical (Lefloch V Lazareff 1994; Kessel-Deynet V 


Burkert 2003 Miao et al.||2006[ |2009 ). Observations have 


also reported several type B/C BRCs with wide tail mor¬ 


phological structures, e.g. SF074 ( [Thompson et al.| |2004 
Kusune et al.||2014 ). 


In summary, a variety of non-standard morphological 
structures can be obtained by changing the inclination of 
cloud G1200(5). At low inclination, linear fragment-core 
structures form, at middle inclinations, asymmetric type B 
and C structures, and approaching 90°inclinations these re¬ 
turn to symmetrical type B and C structures. 


3.2.2 The location of high density cores 

Figure [7] illustrates the distribution of the x-displacements 
of condensed cores formed in the 7 simulations of different 
initial inclination angle p. A point to note is that, where 
cores form in very close proximity, their positions cannot be 
distinguished on the scale used, regardless of the size choice 
of the marker. It can be seen that when p > 0°, core(s) 
always form at the side of positive x axis and the change of 
the x displacement with p follows the trend of x p with <p. 
This illustrates an extremely high preference for high density 
core(s) to form around the apex x p in each simulation. 



(p, inclination 


Figure 7. The x displacement of the cores formed in G1200(5) 
with different initial inclination angles. The short indicates 
Xmax, ‘+’ x p , ‘o’ for the cores with 10 6 ^ n pea k ^ 10 12 cm -3 , 
and V for the extremely high density cores with n pea k > 10 12 
cm -3 . 



Time (Myr) 


Figure 8. The evolution of maximum density of the core(s) in 
the cloud G1200(5) with different initial inclination angles. 


3.2.3 RDI triggered core formation efficiency 

In order to examine the effect of varying inclination angle on 
the efficiency of RDI triggered high density core formation, 
we investigate the angle dependence of the core formation 
time and the accumulated core mass. Figure [8] shows the evo¬ 
lution of the maximum density in the 7 simulations with dif¬ 
ferent inclination angles. In general, the evolution of n max in 
each simulation passes through three different phases: quasi- 
linear, quasi-stable and steeply rising, although the time of 
each phase varies with inclination angle p. These phases may 
correspond to the different stages of RDI triggered high den¬ 
sity core formation. Beginning with initial compression by 
an RDI induced shock, mass accumulates to form the con¬ 
densed region, which finally collapses to form high density 
core(s). The time for a simulation to reach the highest den¬ 
sity (high density core formation) decreases with inclination 
angle, from 0.19 down to 0.086 Myr for simulations with 
p — 0 and 90° respectively. This is because the curvature 
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Figure 6. The final snapshot of the column density from six simulations of G1200(5) with different initial inclination angles; left to right 
then up to bottom are per: 0, 15, 30, 60, 75 90 °. 



Figure 9. The variation of the total core mass of the cloud 
G1200(5) over different initial inclination angles. The symbols 
and their represented properties of cores are shown in the legend 
box. 


of the EUV radiation illuminated surface around the apex 
of the prolate cloud increases with p, this is beneficial for 
mass accumulating towards a geometric focus underneath 
the curved surface. This decreases the path length to the 
focus, and so also the time for accumulating mass around 
it and triggering gravitational instability of the condensed 
region; albeit with less material having been ‘swept up’ by 
the front. Due to this shortening path length, the total mass 
of the high density core(s) is therefore expected to decrease 
with p. 

The relation between the total high density core mass 
and inclination angle p is shown in Figure [9] which reveals 
that the total high density (n ^ 10 6 cm - ) core mass is 
about 4.1 M 0 , in simulation of p = 0, then decreases with 


p, becoming 0.4 M 0 at p — 90°. This is consistent with the 
above expectation. However the total mass of the extremely 
high density core(s) (n > 10 8 cm -3 ) doesn’t show an ob¬ 
vious similar relation with (p to the mass of high density 
core(s). This may be for several reasons. Firstly, as the sim¬ 
ulation is curtailed before the end of all of the dynamic pro¬ 
cesses, extremely high density cores may be on the cusp of 
formation, but have not yet had the opportunity. Secondly, 
a large degree of random variation is expected with such 
gravitational instability-based collapse. Accumulating mass 
in compact high density regions can be fairly consistent, 
but the collapse of these regions or parts of these regions to 
extremely high densities through gravitational instability is 
chaotic and varies strongly with the exact configuration of 
the dynamics. 

From the above result, it is understood that as a cloud 
is rotated from p = 0 to 90°, the morphology of the cloud 
varies from fragment-core filamentary structures to varied 
asymmetrical BRCs, then to standard symmetrical BRCs as 
p approaches 90°. At high p, RDI triggered star formation 
occurs earlier with lower masses at high densities. However, 
for the sampling performed, such pattern is not obvious for 
the extremely high density core formation, which appears to 
be chaotic. Next we look at the effect of the initial shape of 
an inclined prolate cloud on its final morphological structure 
and RDI triggered star formation efficiency. 


3.3 Evolution of G1200 series with 1.0 ^ 7 ^ 8.0 

The initial shape of the clouds investigated here is defined by 
the ratio of the major to minor axis which are in the range of 
1.0 ^ 7 ^ 8.0. When the initial density and total mass are 
constants, the major axis changes with 7 . The inclination 
angle is kept at 45°. The ionization radiation penetration 
depth ranges from 0.225 to 0.320% as listed in Table [5] 
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Ratio, 7 


Figure 11. The x displacement of condensed core in G1200 series 
with an inclination of ip = 45% for varied 7 . The short £ -’ indicates 
Xmax, ‘+’ x p: ‘o’ for the cores with 10 6 ^ %> e ak ^ 10 12 cm -3 , 
and V for the extremely high density cores with n pea k > 10 12 
cm -3 . 


3.3.1 Morphology of final structures 


The dynamical evolution of these clouds is also qualitatively 
similar to that discussed in Section 10 We focus on the 
analysis of the effect of the initial shape of a cloud on the 
final morphology, the location(s) of the condensed core(s) 
and the total core masses. Figure [To] presents six simulation 
result of cloud G 1200 with different initial major to minor 
axis ratio 7 =1.0, 1.5, 2.0, 2.5, 3.0, 5.0, selected from simula¬ 
tion results for 19 different initial clouds, described in Table 
[2] With increasing 7 , the final morphology changes from a 
standard (symmetrical) type A (7 = 1.0), to an asymmetric 
type A BRC (7 = 1.5), then to an asymmetric type B BRC 
(7 = 2.0), next to evolve to an asymmetrical type C BRC 
(7 = 2.5). With further increase (7 ^ 3.0) only a filamentary 
structure can form. 

In these six simulations, high density cores form at the 
head of the structure except in the case of 7 = 2.5, where the 
highly condensed core is embedded in the small ‘nose’ struc¬ 
ture just outside the head. It is interesting to notice that the 
morphological structure of SF046 (CGI) from observations 
(Harju et al. 1990 Haikala et al.||2010| |Makela &; Haikala 


2013) bears a similar feature to that of 7 = 2.5. We leave the 


discussion for the formation mechanism of this morphology 
to the Section |3.4| where we present and analyse more re¬ 
sults of various kinds of morphological structures. It is also 
emphasised that we don’t intend to compare the detailed 
physical structures between simulation and observation on 
SF046 in this paper, but only to point out the link between 
the interesting features from these simulations to such obser¬ 
vations. Detailed comparison of physical properties requires 
further simulations, and will be addressed in a subsequent 
paper. 


3.3.2 Displacement of the high density cores 


Figure 11 displays the x-displacement of condensed core(s) 


with changing geometry, based on the simulation results of 
19 clouds with different 7 . 

With increasing 7 , x p moves further in the +x direction. 


Figure 12. Total mass of condensed cores in cloud G1200 series 
at an inclination of 45% for varied axial ratio 7 . 


This pattern is displayed in the same form for the condensed 
cores. A very clear trend of the core formation occurring at 
or around the apex for all ratios is exhibited. 


3.3.3 The total core mass 


The trend of total mass of both the high density core(s) 
(empty circles) and extremely high density cores (solid black 
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circles) is a decrease with 7 . This is illustrated in Figure 
The total mass of high density cores has an approximate 
pattern of maximum mass of 2.0 Mq at 7 ~ 2.25 (with an 
outlier as high as almost 3.0 Mq at 7 = 2.75) and a minimum 
of 0.6 Mq at 7 = 8 (with an unusually low value of only 
around 0.2 M 0 for 7 = 6.0). Clouds with lower ellipticity 
appear to have a higher probability of collecting greater mass 
in the final cores. This can be expected as the mass per unit 
length along the major axis decreases with 7 . As a result, 
the amount of mass able to collapse toward the apex region 


decreases, which is similar to the case when ip — 0° (Kinnear 
et al.|2014 ). When 7 ? ^ 0°, due to the wide range of variation 
in the final morphologies of the clouds, the amplitude of the 
fluctuation is similarly increased. 

From the same figure, we can see the mass of extremely 
condensed cores (solid black circles) in a cloud also decreases 
with 7 in a similar way as the high density core. The mass 
of these cores corresponds to the precursor mass for possi¬ 
ble proto-stars forming later in the evolution of the cloud, 
having a maximum of ~ 1.6 and a minimum of ~ 0.2 Mq. 


3.4 Morphological evolution of clouds of different 
initial density and under different ionizing 
flux 

After the investigation of the effects of both initial shape 
and inclination angle of a prolate cloud on its RDI triggered 
dynamical evolution, we now explore the consequences of 
changing the initial density of a cloud and the ionizing ra¬ 
diation flux in the morphological development of a cloud of 
mass 30 Mq. This is essentially to examine the consistency 
of the ^euv parameter in characterising a wider variety of 
cloud parameters. 

A sequence of 8 different initial densities are chosen as 
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Figure 10. The final snapshots of the column density from six simulations of G1200 with 7 equal to the values shown in the top-left 
corner of each panel. The inclination angle is 45° for all of the initial clouds. 


Index 

n 

a 

4uv(^i) 

^EUv(^2) 


cm -3 

pc 

% 

% 

G100 

100 

1.797 

11.779 

23.558 

G200 

200 

1.427 

3.708 

7.416 

G400 

400 

1.132 

1.17 

2.34 

G600 

600 

0.989 

0.594 

1.188 

G700 

700 

0.940 

0.460 

0.920 

G800 

800 

0.899 

0.368 

0.736 

G1000 

1000 

0.834 

0.254 

0.508 

G1200 

1200 

0.784 

0.187 

0.374 


Table 4. Parameters of molecular clouds of mass 30 M®, an 
inclination angle cp = 60°, 7 = 2 , and a cr it = 0.079 (pc). From 
left to right, columns are the cloud index, initial number density, 
the major axis and d EE v defined by Equation [ 8 ] with different 
ionization flux F\ and F 2 being 1.0 X 10 9 and 2.0 X 10 9 cm -2 s -1 
respectively. 


100, 200, 400, 600, 700, 800, 1000, and 1200 cm -3 . Two dif¬ 
ferent ionization fluxes of F\ = 1.0 x 10 9 and F 2 = 2.0 x 10 9 
cm -2 s _1 are applied for a total of 16 simulations. An incli¬ 
nation angle of 60° and 7 = 2.0 is used for all cases. This 
is as a result of the simulation of G 1200 (5) exhibiting the 
formation of the ‘nose’ structure in the vicinity of these ge¬ 
ometries. As such, we can also expect to get some clue of how 
the ‘nose’ structure is developed and changes with the initial 
properties of a cloud in this set of simulations. The corre¬ 
sponding parameters of the clouds are listed in Table [4] It 
is seen that the change in EUV radiation penetration depth 
over the initial density (100 ^ n ^ 1200 cm -3 ) is almost two 
orders of magnitude, a much more dramatic variance than 
from changing 7 or (^. 


3.4-1 An overview of the final morphological structures 

Figure [13] displays the final snapshots of the column density 
of the 8 simulations with ionization flux of F\. The mor¬ 
phology of the final structure of a cloud changes from an 
asymmetrical type C BRC to a filamentary structure and 
then to an irregular structure, as the initial density is de¬ 
creased from 1200 to 100 cm -3 , which leads to an increase 
of d E uv from 0.187 to 11.779%. 

As shown in the top row panels, an asymmetrical type 
C BRC is developed in G 1200 , G 1000 , G800, and G700, 
with d E uv of 0.187, 0.254, 0.368 and 0.460% respectively. 
Although they all bear an asymmetric type C BRC mor¬ 
phology, the minor axis becomes narrower with increasing 
cfeuv, which results in more cloud material being photo- 
evaporated from their star facing surfaces. When the initial 
density n — 600 cm -3 (cfeuv = 0.594%), the prolate cloud 
evolves into a filamentary structure, as shown in the first 
panel in the bottom row of Figure [T3] Clouds of initial den¬ 
sities of 400, 200 and 100 cm -3 (d E uv = 1.17, 3.708 and 
11.779 % respectively) are all seen to form irregular struc¬ 
tures. 

The final morphological structures from the simulations 
using doubled ionization flux F 2 are presented in Figure [14] 
Similar morphological structure sequences to those seen in 
Figure |T3] are observed. An asymmetrical type C BRC is 
formed from G 1200 , G 1000 and G800, a filamentary struc¬ 
ture from G700, and irregular structures in G600, G400, 
G200 and G100. The only difference is that the morphologi¬ 
cal transition point shifts to a cloud of higher initial density 
compared to the F\ simulations. 

The first transition from type C BRC to filamentary 
morphology occurs around G600(Fi) and G700(F 2 ), and the 
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Figure 13. The final snapshot of the column density of clouds of initial mass of 30 Mq 7 = 2.0 and tp = 60° but different initial 
as shown in the top-left corner of each panel. The ionization flux — 10 9 cm -2 s —1 
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Figure 14. The final snapshot of the column density of the same 8 clouds as in Figure [13| (initial mass of 30 Mq 7 = 2.0 and cp = 60° 
but with increased incident flux of 2 x 10 9 cm -2 s —1 . 
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Figure 15. The core formation time of the F\ and F 2 group 
clouds vs their initial densities. 



Initial density (cm 3 ) 

Figure 16. The total mass of the cores formed in F\ and F 2 
group clouds vs their initial densities. 

second transition from filamentary to irregular morphology 
is around G400(Fi) and G600(F2). These appear to corre¬ 
spond to similar cfeuv in each case, for both fluxes. The 
former transition occurs at cfeuv ~ 0.8 and the latter at 
cfeuv ~ 1.1. This indicates that cfeuv remains a characteris¬ 
tic parameter for the dynamical evolution of a cloud for the 
parameter space so far investigated. 

Figures [15] and [16] present the core formation time and 
the total mass of the RDI triggered condensed cores vs the 
initial density of the clouds of the above two groups of simu¬ 
lated clouds, with black lines for Fi group and red for F 2 . It 
is shown that the core formation time in both group clouds 
decreases with the initial density of a cloud. This is because 
of the stronger initial gravitational bonding in the high den¬ 
sity cloud. For the same reason, the total mass of the con¬ 
densed cores follow an increasing trend with the initial den¬ 
sity of the cloud. However the effect of the strength of the 
EUV radiation flux on the efficiency of triggered core forma¬ 
tion vs initial density is two fold. When the initial density 
is below 600 cm -3 the clouds illuminated by the higher flux 
(F 2 ) take a longer time to form condensed cores, and col¬ 
lect slightly more mass in their condensed cores than in the 
clouds illuminated by the lower flux (Fi). This trend is re¬ 


versed when the initial density is higher then 600 cm -3 , i.e., 
clouds illuminated by higher flux take less time to form con¬ 
densed core and collect less mass into their condensed cores 
than the clouds illuminated by lower flux. 

Overall, the information delivered in the above two 
figures suggests that the RDI triggered star formation ef¬ 
ficiency is much higher (shorter core formation time and 
higher total core mass) in molecular clouds of shorter cfeuv 
(corresponding to higher initial density) in both Fi and F 2 
group clouds. 

From the above two sets of simulations, we can see that 
when (few ^ 1.0%, the deep penetration of the EUV ra¬ 
diation into the cloud causes irregular structure formation. 
A closer look at the above two morphological variation se¬ 
quences tells that with increasing (few, the tiny ‘nose’ struc¬ 
ture around the apex of the asymmetrical type C BRC grad¬ 
ually grows and becomes an irregular ‘horse-head’ structure. 
In order to understand how the ‘nose’ and ‘horse-head’ struc¬ 
tures form, we further investigate the evolutionary process 
of one of the above clouds, which presents a clear horse- 
head-like morphology. 


3.4-2 Growing from ‘nose’to ‘horse-head’ structure 

We chose the cloud G400(F2) to analyse its evolutionary pro¬ 
cess due to its formation of a particularly distinctive ‘horse- 
head’ structure. Figure p~7| describes a six-snapshot column 
density evolutionary sequence of cloud G400(F2), which 
shows a growing process from a small ‘nose’ to a ‘horse- 
head’ morphological structure. The corresponding RDI in¬ 
duced shock velocity profiles for each of the snapshots are 
plotted in Figure [18] which reveals the kinematic mechanism 
for the morphology formation in G400(F2). 

From Figure [TT| we can see that a shocked layer forms 
at t — 0.046 Myr, which surrounds most of the surface of the 
front semi-ellipsoid and a small part of the top surface of the 
back semi-ellipsoid. The maximum magnitude of the shock 
velocity v s = 12 km s -1 , calculated as the peak velocity of 
non-ionized gas within the condensed shock front. Most of 
the shocked layer surrounding the front semi-ellipsoid prop¬ 
agates into the cloud smoothly toward (+£, -z) direction. 
Very close to the apex, the shock velocity of the front semi¬ 
ellipsoid is toward the direction of (x, —z), and that of the 
back semi-ellipsoid is (—x,—z), the gas material from the 
two sides of the apex converge beneath the apex to make 
a small ‘nose’ structure at t — 0.084 Myr, as shown in the 
top-right panels Figures [Tt] and [TS] 

While a small ‘nose’ structure forms at the point of (x, 
z) ~ (0.2, 0.2) pc and with v x « 0, below the nose structure, 
the RDI induced shock keeps pushing the surface layer of the 
front semi-ellipsoid moves toward the (+£, —z) direction at 
the shock velocity. This makes the below-nose region of the 
BRC gradually produce a displacement along x direction 
relative to the nose. The nose becomes a short horizontal 
fragment at t = 0.13 Myr, in which the direction of the RDI 
induced shock velocity is mainly along z direction. The short 
horizontal fragment gradually forms a ‘hook’ at t — 0.17 Myr 
as shown in the lower middle panels of Figures [T7| and [18] 
Further propagation of the shocked layer below the hook 
toward (+x, —z) direction makes the ‘hook’ structure grow 
to look like a ‘horse-head’ as shown in the last panels in 
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the above two figures. The final ‘horse-head’ structure bears 
some resemblance to the well known horse-head nebula. 

It is seen from the above analysis that the appearance 
of a ‘horse-head’ morphology in this instance is a conse¬ 
quence of growing a small ‘nose’ structure. The ‘nose’ struc¬ 
ture forms if the gas around the apex of an ellipsoidal cloud 
converge beneath that apex due to a strong RDI induced 
shock effect, as described above. This requires that the cloud 
should have a reasonably high inclination angle, i.e., tp > 45° 
in our simulations. However the formation of a ‘nose’ struc¬ 
ture does not necessarily lead to a ‘horse-head’ structure. 
For example, a small ‘nose’ formation is seen from the fi¬ 
nal morphologies of G700-1200 in the F\ series and G800- 
1200 in the F 2 series. Whether the ‘nose’ can grow into a 
‘horse-head’ also depends on the initial ionization penetra¬ 
tion depth, i.e., dsuv > 1.0%. The latter is to guara ntee 


a higher enough shock velocity ^ 'v s ~ FL ^ Tian ~^ (Bertoldi 
1989), which can push the majority of the gas below tEe 


nose toward the (+x, —z) direction. 


3.5 Development of other irregular morphologies 

Further simulations were conducted, to explore if there are 
other structures which form different irregular morphologies 
to a ‘horse-head’ as a consequence of the variation of initial 
conditions of the cloud and ionization flux. A series of 16 
simulations with G400(F2) of different 7 = 2.0, 2.5, 3.0, 3.5 
and tp 8=3 60,70,80,85° are conducted. For all 16 clouds, 
^euv > 1.0%. In this section, we concentrate on the de¬ 
scription of the morphological variations, as the underlying 
mechanisms appear to be similar, and they are of interest 
due to their variety rather than exact evolutions. 

The first column panels in Figure [T9] present the final 
morphologies of G400(F 2 ) of 7 = 2 and 4 different initial 
inclination angles, which show that when the initial angle 
tp increases from 60°, the ‘horse-head’ in the final structure 
becomes longer in z direction. The connection point between 
the ‘horse-head’ and the below-horse-head part moves down¬ 
ward with (p. At 85 , the whole morphology is an ‘elephant 
trunk’ like structure. This is because the ‘nose’ structure 
formed in the first tens of Kyr of the evolution increasingly 
stretches along the z axis with increasing ip. The first row 
panels in Figure [19] describe the morphology change of cloud 
G400(F2) of inclination angle tp = 60° and varied 7 . When 
the initial prolate cloud becomes more elliptical, in addition 
to a ‘horse-head’ structure, the RDI induced shock also trig¬ 
gers linear gravitational instability along the direction of the 
major axis in some clouds. 

The 16 panels in Figure [19] present 16 different final 
morphologies which appear in the clouds of the initial ( 7 , tp) 
as stated above. We can conclude that the larger the incli¬ 
nation angle, the longer the ‘horse-head’ sub-structure; and 
that the more elliptical the cloud initially, the thinner the fi¬ 
nal morphology is. Although, as in previous instances, their 
structural details differ across the range of 7 . 

Through viewing more simulation results obtained us¬ 
ing different combination of (Flyman, n, 7 , tp), we find that 
the variety of the final morphology is extensive for ^euv > 
1 . 0 %, although it is impossible to show all of them in this 
paper. 


4 CONCLUSIONS 

The results from the simulations conducted show that uni¬ 
form density prolate clouds at an inclination to a plane par¬ 
allel EUV radiation field can evolve to a variety of morpho¬ 
logical structures often found at various Hu boundaries. The 
parameter of EUV radiation penetration depth, cZeuv, of a 
given cloud can be a good indicator as to whether it would 
evolve to a regular BRC or irregular structures. 

The final morphology of the clouds of initial denisty 
(1200 cm -3 ) and ellipticity (7 = 2), and with efeuv < 1.0% 
changes from a filament to an asymmetrical BRC, then to a 
symmetrical BRC, as the inclination angle tp increases from 
0 to 90°. The core formation time and the total mass of 
the RDI triggered condensed cores both decrease with tp. 
Therefore EUV radiation triggered star formation is quicker 
in high tp prolate clouds, but the total mass of the triggered 
stars is higher for low tp cloud. 

For clouds of fixed initial density (1200 cm -3 ) and in¬ 
clination angle (tp = 45°), as 7 increases from 1 to 8 , the 
final morphologies changes from a symmetrical to asymmet¬ 
rical BRC, then to filamentary structure. RDI triggered star 
formation efficiency decreases with the increase of 7 , due to 
a greater quantity of gas material being evaporated with a 
high cross sectional area to the incident radiation. 

It is found that the final morphology of a prolate cloud 
is very sensitive to the initial density. For fixed initial shape 
(7 = 2 . 0 ) and fixed inclination angle (tp = 60°), the mor¬ 
phology of a cloud changes from an asymmetric type C BRCs 
to a filament then to irregular structures, such as a ‘horse- 
head’, as the initial density decreases from 1200 to 100 cm -3 . 
EUV radiation triggered star formation efficiency is found 
to be higher in regular BRCs (with cfeuv < 1.0%) than in 
irregular structures (with cfeuv > 1 . 0 %). 

Based on the above systematic investigation on the evo¬ 
lution of molecular clouds (or clumps) at an Hu bound¬ 
ary, we can now suggest a unified formation mechanism for 
the different morphological structures seen in the published 
images of a variety of Hu regions: Giant molecular clouds 
(GMCs) are clumpy and the numerous clumps are of dif¬ 
ferent shapes, from spherical to highly elliptical, and with 
their major axes not exhibiting a strong preference for any 
particular alignment within their environment. When a star 
or a cluster of stars forms inside a GMC, the interaction be¬ 
tween the ionizing radiation with the surrounding molecular 
clumps of different shapes and orientations creates struc¬ 
tures of various morphologies. The spherical and prolate 
clumps aligned with the radiation flux form standard type 
BRCs, the inclined prolate clumps with a shallow ionizing 
radiation penetration depth form asymmetrical BRCs or an 
inclined linear structures; the inclined prolate clumps with a 
deep radiation penetration depth form a variety of irregular 
structures, the well known ‘horse-head’ structure is only one 
of them. 

In our next paper, we are going to investigate the evo¬ 
lution of physical properties of a few observed asymmetrical 
BRCs and the RDI triggered star formation inside them. 
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Figure 17. The column density evolution of a cloud G400(-F2), with an inclination angle of 60°. 
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Figure 18. The evolution of the RDI induced shock velocity of cloud G400(F2) overlapped in the corresponding cross-section density 
snapshot. The magnitude of the velocity is between 2 and 12 km s —1 . 
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Figure 19. The final morphologies of the clouds G400(T2) of different initial geometry and inclination angle. The axis scales vary 
between panels, because of the different initial scale of the major- and minor axes for the clouds depicted. 
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